Electronic properties of graphene nanoribbons under gate electric fields 
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Quantum-dot states in graphene nanoribbons (GNR) were calculated using density-functional 
theory, considering the effect of the electric field of gate electrodes. The field is parallel to the GNR 
plane and was generated by an inhomogeneous charge sheet placed atop the ribbon. Varying the 
electric field allowed to observe the development of the GNR states and the formation of localized, 
quantum- dot- like states in the band gap. The calculation has been performed for armchair GNRs 
and for armchair ribbons with a zigzag section. For the armchair GNR a static dielectric constant 
of £ ~ 4 could be determined. 
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I. INTRODUCTION 

After their experimental realization, single and bilayer 
graphene sheets and nanoribbons have attracted inten- 
sive attention due to their peculiar properties, which 
make graphene and its derivatives one of the mostpromi- 
nent material classes for future nanoelectronics.^^ The 
vast range of^ossible applications is due to the high car- 
rier mobility,^ ^ and remarkably long spin lifetimes and 
phase coherence lengths, which are particularly valuable 
for quantum information processing.^ The peculiar elec- 
tronic structure of gapless semiconductor graphene pre- 
vents electrostatic confinement due to Klein tunneling 
and gives rise to unique transport properties.^ Experi- 
mentally, graphene-based tunable nanodevices can be re- 
alized, as shown, e.g., for graphene nanoribbons ,'I^H^in- 
terference devices and graphene quantum dots.^^^^ 
Graphene nanoribbons (GNRs) play a particularly im- 
portant role since they are often used in field-effect tran- 
sistors (FET) setups where the current flow is regulated 
by a gate electric field. Additionally, gated nanorib- 
bons can also be used to create quantum dots. Hereby, 
the spin qubits in quantum dots are a promising can- 
didate for quantum information processing, for which 
graphene is better suited than III-IV semiconductors 
due to its reduced hyperfine coupling and spin-orbit 
interaction.'^^^ For nanoelectronics applications, the 
preferred width of a GNR is in the range of 1- 3 nm as 
for wider ribbons the band gap diminishes. '^^'^ 

For a proper theoretical description of GNR FET and 
GNR-based quantum dots, both the effect of an exter- 
nal electric field and the precise electronic structures of 
the ribbons have to be considered. In ribbons of small 
width the details of the termination and exact structure 
of the edges play an important role for ribbon's elec- 
tronic structure due to the quantum confinement effects, 
which have to be properly taken into account. The two 
most common types of edge termination of the ribbons 
lead to very different electronic structures in the vicinity 
of the Fermi energy: while armchair ribbons are insu- 



lating with a band gap of up to approxim ately 2.5 eV, 
zigzag ribbons exhibit metallic edge states.'^^'^ For the 
former, the band gap decreases with increasing width of 
the ribbon and oscillates in magnitude with a periodicity 
of three as more rows of carbon atoms are added.E21 Ex- 
perimentally, armchair ribbons with the width of 1 nm 
with well-defined edges^^ and FET with widths down to 
2 nm have been realized. On the other hand, ribbons 
with the width of around 50 nm have been used to create 
quantum dots.^^ 

The sensitivity of the electronic properties of narrow 
GNRs to an external electric field and to the details of the 
structure, as well as the mechanism for the formation of 
the quantum dot states, calls for a description and anal- 
ysis within the highly accurate density-functional theory 
techniques (DFT). In this work, we study the properties 
of narrow GNRs in an electric field, generated by finite 
gates, by employing a new scheme specifically designed 
for this purpose and realized within the full-potential lin- 
earized augmented plane-wave (FLAPW) method. In 
contrast to common approaches, the latter scheme allows 
to consider any given distribution of charge in the gates, 
thus enabling the versatility necessary for studies of com- 
plex nanostructured devices from first principles. Here, 
we describe the details of this approach as implemented 
within the FLAPW method, and study the develop- 
ment of the electronic structure and screening properties 
of narrow GNRs as the strength of the electric field is 
varied. We focus particularly on the emergence of the 
quantum-dot like states, localized under the gates, and 
analyze their spatial distribution and symmetry proper- 
ties. 



II. INCLUSION OF THE ELECTRIC FIELD 

In a realistic device the electronic properties of the 
graphene nanoribbon will be modified by gate electrodes. 
While more complex interactions between the gates and 
the graphene nanoribbon exist, the most prominent ef- 
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feet of the gates is the electric field they generate. These 
electric-field effects can be included in a calculation based 
on DFTpIH^ and in the past several different imple- 
mentations of external fields in DFT codes have been 
reported I^^JESI However, in contrast to most of these ap- 
proaches, which model a situation in which the electric 
field is applied normally to a surface, the key effect of the 
gates we study is the field and potential distribution in 
the plane of the graphene ribbon. 



A. Computational scheme 

We use the film FLAPW method^^ as implemented in 
the FLEUR code^^ for our calculations. Two slightly dif- 
ferent schemes to include the electric fields of the gates 
have been tested in which the effect of the gates is mod- 
eled by: (i) a charged sheet far in the vacuum in which 
the charge varies parallel to the GNR or (ii) a plane in 
the vacuum at which a varying boundary condition to the 
Coulomb potential is applied. While the second approach 
is closer to the picture of a metallic gate applying an elec- 
tric field, the generation of the potential from a charge 
distribution is performed in all existing DFT codes and, 
hence, the first approach is easy to implement. We dis- 
cuss all details of the implementation of the gate electric 
field in the appendix. 

For modeling the effect of a gate electrode on a 
GNR we have chosen a free-standing 13-carbon- atom- 
wide graphene nanoribbon of armchair type [see Fig. 
[TJa)]. The dangling bonds are passivated using hydro- 
gen leading to a ribbon which is about 5.1 nm long and 
1.7 nm wide, consisting of 312 carbon and 48 hydrogen 
atoms. The setup is periodically repeated in x-direction 
to form an infinite ribbon. As our code requires a two- 
dimensional periodicity, we also have to repeat the ribbon 
in y direction. A supercell approach with a separation of 
5 A between adjacent cells is used to simulate a single 
isolated GNR. The center of the GNR is sandwiched be- 
tween positively charged top and back gates whose charge 
is compensated by negatively charged top and back gates 
such that both the GNR and the two charge sheets re- 
main as a whole charge neutral. The positively charged 
gate has a size of 25.6 x 12.3 and the distance be- 
tween the GNR plane and the charge sheets is 7 A. For 
simplicity, we used a homogeneous surface-charge-density 
distribution on the gates [method (i)]. However, for the 
dimensions of this system, the potential at the position 
of the plane of the GNR is qualitatively the same for a 
gate with constant potential - the potential only differs 
strongly close to the gate (not shown). The calculations 
were performed with the Perdew-Burke-Ernzerhof (PBE) 
exchange-correlation functional.E^ 



B. Screening of the electrical field 

The height dependence {z dependence) of the electric 
field is shown in Fig. [TJb), cutting vertically through the 
GNR (with y being at the middle of the ribbon). In the 
plane of the ribbon (z = 0), the field is fully in plane and 
parallel to the x direction. Above the plane (\z\ > 0), 
an out-of-plane component exists, which grows towards 
the center {x = 2.55 A) and the left /right side (x = 
and X = 5.1 A). At exactly those points, the field is 
perpendicular to the GNR plane and the field strength 
is low, while for the points in between the field is nearly 
in-plane and stronger. Thus, the potential should cause 
a charge accumulation under the gate. 

Figure [2|a) shows the gate potential along the center 
of the ribbon for a gate charge of 0.0016e/nm^. These 
curves have been obtained by subtracting the total po- 
tential without an applied field from the potential with 
applied gate field canceling the potential from the mostly 
inert ion cores. The (in-plane) applied electric field at the 
position of the ribbon is given by the slope of the solid 
black curve and has a maximal value of 0.056 eV/A. The 
initial potential due to the applied electric field (black 
solid curve) gets reduced (red dashed curve) due to the 
screening by the electrons in the ribbon. This screen- 
ing effect is determined by evaluating the potential of 
the self-consistent calculation. From the reduction - the 
blue dotted curve has been obtained by multiplying the 
red dashed curve by 3.8 to match the black solid curve 
- a static {uj = 0) dielectric constant £ ~ 3.8 can be de- 
duced for the center of the ribbon. As the calculation has 
periodic boundary conditions (in x-direction) this static 
dielectric constant £{q) corresponds to a wavevector of 

qo = 0.123 A . The value of e{qo^uj = 0) ^ 3.8 agrees 
with the GW results of Ref. SSI One should note, that 
the static dielectric constant has been evaluated in the 
linear response limit of a small electric field, for stronger 
fields and especially close to the crossover of states [cf. 
Fig. [4|b)], different values are obtained. 

The dependence of the screening on the distance from 
the GNR plane is shown in Fig. [2|b); £ reduces only 
slightly from 3.8 to 3.4 over 1 A. This relatively constant 
screening corresponds to presence of change density due 
to the Pz orbitals, which is still significant at these dis- 
tances. For larger distances to the GNR e decreases ex- 
ponentially, approaching the vacuum value of 1. 



III. RESULTS 

A. Straight armchair ribbon 

In this section, we focus on the formation of quantum 
dot states in an armchair GNR with a width of 13 carbon 
atoms by applying an electric field as outlined in the 
previous section. The edges of the ribbon are hydrogen- 
terminated and it is placed between the electrodes as 
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FIG. 1. (Color online) (a) Free-standing, hydrogen- 

terminated armchair graphene nanoribbon of 13 atoms width. 
The blue rectangles indicate the location of the positively 
charged top and bottom gate electrodes. Their charge is com- 
pensated by the negatively charged top and bottom electrodes 
(shown in red), (b) Side view showing the electric field cre- 
ated by the electrodes (for y at the middle of the ribbon). 



visualized in Fig.[TJa). If one compares the total electron 
density of this structure without field and with a field of 
0.67 V/A (this corresponds to a charge of ±0.19e/nm^ on 
the gates), the charge accumulation and depletion can be 
visualized by a density-difference plot as shown in Fig. |3] 
Predominately charge of the orbitals is localized under 
the positively charged gate forming a "quantum dot" of 
about 20 A diameter. 

To investigate the effect of the electric field on the 
electronic structure of the GNR, we first consider in 
Fig. [4]^ a) the bandstructure of the ribbon without field. 
The bandgap is about 0.9 eV and bounded by valence 
and conduction band states at the F point. Introduc- 
ing an additional periodicity in x-direction by the field 
leads to a backfolding of the bandstructure along F-X. 
If the field is modulated with a periodicity of 51 A (12 
unit cells in x-direction), the backfolding occurs at the 
red lines shown in Fig. p[a). The Brillouin zone extends 
then from F to X^ In the following we focus on the 
states at the F point, since they form the valence and 
conduction band edge. 

Switching on the electric field leads to a splitting or a 
shift of the eigenenergies as shown in Fig.[4]^b): The two- 
fold degeneracy of the states at —0.75 eV and +0.7 eV, 
which results from the backfolding of the Brillouin zone, 
is lifted and the bands split almost symmetrically, both 
for negative and positive fields. The slight asymme- 
try between positive and negative electric field is intro- 
duced by the shape of the gate electrode as shown in 
Fig. [ija). If there were just a one-dimensionally mod- 
ulated field acting on the GNR, this asymmetry would 
vanish. The splitting can then be understood from 
the consideration of a system with just two eigenstates 




FIG. 2. (Color online) (a) Potential due to the gate elec- 
tric field along center of the GNR {z = 0). The black solid 
curve shows the applied gate potential, the red dashed curve 
the screened potential. The blue dotted curve is obtained by 
multiplying the red dashed curve by 3.8, indicating a static 
dielectric constant of s ~ 3.8. (b) Dependence of the static 
screening on the distance from the plane. 




FIG. 3. (Color online) Total charge density difference be- 
tween a calculation with gate electric field and without; the 
density is reduced (blue) at outside and increased (red) in the 
center. 



with energies Cq and Cq , perturbed by a periodic po- 
tential with Fourier coefficients Vi. At the i^^ backfold- 
ing of the bandstructure, the states split according to 

eiAV) = ^(4 + eo ) ± ^JI{4-^o ? + \V^\^■ As the 
states are degenerate, the splitting is linear in In 
case of interacting states that are non-degenerate, the 
evolution of these states starts parabolic at small fields 
and evolves into a more linear behavior at stronger fields 
(T/i»(e+-eo-)). 

With increasing field strength the uppermost valence 
band states move up in energy, towards the Fermi level, 
while the lowest conduction band states move down by 
a similar amount. Their interaction causes again a de- 
viation from the linear behavior and at a field strength 
corresponding to 0.4 V/A the first crossover of states oc- 
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FIG. 4. (Color online) (a) Band structure for the smallest, 
13 carbon atom wide GNR (Y12 as long as the gated GNR) 
with four staggered carbon atom in x direction and totally 
26 carbon atoms; the zero- fie Id band gap of about 0.9 eV 
is shaded in yellow, (b) Single-particle eigenenergies gated 
GNR for the F point in dependence of the applied electric 
field. The energy is relative to the Fermi energy; the abscissa 
shows maximal electric field (gradient of the potential before 
screening) at position of the ribbon. The states at the dashed 
vertical line [(a)-(h)] are depicted in Fig. [s] 




curs, i.e. a valence band state becomes unoccupied while 
a conduction band state is populated. In order to get 
a better understanding of this scenario, we display the 
single-particle states in this energy region in Fig. [5] for 
a field of 0.67 eV/A [vertical line in Figjij^b)]. The first 
crossover occurs for states shown in Fig.|5je) and (h), i.e. 
with state (h) the first quantum dot state gets localized 
under the gate electrode. State (e) has more density out- 
side the positively gated region; it raises in energy with 
increasing field and gets unoccupied. 

At stronger fields (0.72 eV/A) the state Fig. [5|f) gets 
occupied, which can be regarded as the second quan- 
tum dot state. As for atoms, one expects also for quan- 
tum dots states an increasing number of nodal planes for 
higher states with higher energy, which is indeed what 
one observes: The color of the isosurfaces in Fig. |5] de- 
notes the sign of the wavefunction, which allows to locate 
the nodal planes. The first quantum-dot state (h) has at 
the center one vertical nodal plane (dashed green line) 
and is not nodeless as one might have expected. For the 
next localized state (f) one sees in the charge density a 
horizontal nodal plane, while the state (c) has two: a hor- 
izontal and a vertical nodal plane. The state (b) seems 
to have two horizontal and the state (a) one horizontal 
and two vertical nodal planes. 

We have thus shown how under a gate electric field, 
localized states appear in the band gap; those show a 




FIG. 5. (Color online) Charge density of the single-particle 
levels at a field of 0.67 eV/A, ordered descending in energy; 
(a)-(f) are above, (g) and (h) below the Fermi energy. The 
color of the isosurfaces denotes the sign of the wavefunction. 



nodal structure as one expects for a quantum dot. The 
states still show the structure of the underlying lattice 
and are predominately formed by Pz states as one can 
see in the charge density. 



B. Z-shaped ribbon 

Besides armchair nanoribbons, also zigzag ribbons 
exist, which are known for their (conducting) edge 
states. Those are located just above and below the 
Fermi energy. The combination of both types of edge 
states could be useful for future nanodevices; but miss- 
ing atoms or sections with different edge terminations 
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FIG. 6. (Color online) Armchair nanoribbon with zigzag 
section, (a) Edge state at the Fermi energy (no field), (b) 
State above 0.87 eV above the Fermi energy (no field), (c): 
Same state with an applied gate field of 1 eV/A, now located 
at the Fermi energy. 



might also occur involuntarily during the preparation of 
nanoribbons. We have thus simulated an armchair rib- 
bon of width 13 (as above) but with a zigzag section in 
the middle. 

For the calculation, a supercell with length of 51.3 A 
and width of 25.8 A has been used; the positive gate 
charge has been placed at the center over and below the 
zigzag region (16 x 18.8 A^). Figure [g]^ a) shows (for zero 
field) the existence of an edge state in the zigzag sec- 
tion; the highest occupied (shown) and the lowest unoc- 
cupied state have the same density but locally different 
phases. If one now turns on an electric field, the edge 
states gain density in the middle of the ribbon. Addi- 
tionally, states with more density at the outside raise in 
energy, while more localized states get more localized and 
become lower in energy. This is exemplified in Fig. [6] for 
an electric field of 1 eV/A (maximal potential gradient 
before screening; gate charge is 0.026e/nm^): (b) shows 
a high lying state at zero field (0.87 eV above the Fermi 
energy) , which moves down in energy to the Fermi energy 
in the electric field; at the same time its density becomes 
more localized at the zigzag region, i.e. under the gate. 



IV. CONCLUSION 

A scheme to include arbitrary electrical fields has been 
introduced and its implementation for the film FLAPW 
method outlined. The scheme allows for Neumann and 
Dirichlet boundary conditions and for differently shaped 
gate electrodes. Future use could encompass the calcula- 
tion of the effects of an electrical field on adatom on films 
and a wide range other electric-field related phenomena. 

The gate electric field was then applied to an armchair 
graphene nanoribbon; the mostly in-plane field caused a 
charge accumulation under the gate. For small fields, the 
linear response of the electrons in the ribbon to the elec- 
tric field allowed to determine a static dielectric constant 
of £ = 3.8. The field also drove states into the zero- 
field gap; as has been shown, the previously unoccupied 
states entering the gap region are localized under the 
gate. Those states, mostly formed by Pz orbitals, show 
the structure of the underlying lattice. However, they 
also feature a nodal structure as one would expected for 
quantum dots. The states have been obtained from an 
all-electron density-functional theory calculation, which 
makes a comparison to less precise techniques such as 
tight binding interesting; those techniques have the ad- 
vantage that much larger systems can be treated. Addi- 
tionally, the influence of the electric field on GNR consist- 
ing of zigzag with metallic edge states and semiconduct- 
ing armchair sections has been shown, where the electric 
field moves states towards the Fermi energy, which are 
localized under the gate. 
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Appendix A: Electric field in FLAPW 

For periodic systems, the full-potential linearized aug- 
mented plane wave (FLAPW) basis, used for the calcu- 
lations pre sented above, gives an accurate, all-electron 
description-I^SliMIIl Besides the more common FLAPW 
basis with three-dimensional periodicity, a film version 
of FLAPW exists. It has only two-dimensional periodic- 
ity and semi- infinite vacua in the third direction (chosen 
to be 2;)p5I3SIiniIII13 Hereby, the space is separated into 
three regions: the muffin-tin spheres around each atom, 
the interstitial region between the spheres and a vacuum 
region in which the density decays exponentially to zero; 
see Fig. [7[ The basis set used for expansion of the wave- 
function consists of the functions of the following form: 
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gi(G|| +k|| )v^iG^z^ interstitial 
^G||G^(k||,r) = <^ [aG||,G^(k||)iiG||(k||,z) + 6G||,G^(k||)^G||(k||,z)]e'(^ii+^ii)'' vacuum , (Al) 

[a^j^ {y)ui{r) + 6^^(k)ix^(r)]yL(r) /i-th muffin tin 



where G = Gy +(0, 0, G^) is the reciprocal lattice vector, 
the ul are radial functions with as their energy deriva- 
tives, Yl are spherical harmonics, r = r/|r|, and a and 
h are coefficients chosen such that the basis function is 
continuously differentiable across the interstitial-vacuum 
and interstitial-muffin-tin boundary. The || denotes the 
in-plane {x-y) and ± the out-of-plane {z) component. 
While the interstitial extends to \z\ < zi = I)vac/2, the 
perpendicular wavevectors G±_ are defined with regard to 
D > Dvac, i-e. G± = 27rn/Z), to allow for more varia- 
tional freedom. 



1. Neumann boundary conditions 

The idea is to obtain the potential by integrating the 
surface charge density a from infinity, i.e. 




0- 



Muffm tin 
Interstitial 
Vacuum 
Charge sheet 



FIG. 7. (Color online) Geometric setup in the FLAPW 
method with mufiin-tin spheres, interstitial, and two semi- 
infinite vacuum regions in z direction. \z\ = zi = Dvac/2 
marks the interstitial-vacuum boundaries, \za \ is the position 
for the charge sheets. 



V{r) 



£0 J-oo 



(A2) 



with da/dz = p, where p is the charg e de nsity. The 
complication in evaluating the equation ( A2 ) arises from 
the use of different bases in the different regions. The 
electric field S could be included in the calculation via 
its associated potential Vef as an additional term to the 
external potential V^xt; however, in the presented scheme 
it enters as charge density pef- Their relation is given 
by So^'^Vef = SqV ■ E = — Pef, where Eq is the electric 
constant (for cgs replace Sq by l/Air). 

For the vacuum region, the charge p and th e potential 
V can be expanded in a Fourier series, ^ ^^ 1 ^^ * 



Vg,=o{z) = / / p{z^')dz^^dz^ 

^0 Jz Jz' 



- r {z' -z)p{z")dz'. (A5) 

^0 Jz 



For the nonuniform part, Gm 7^ 0, the Green function 



G(^,^') = ^e-l««l-l-^'l 



(A6) 



can be used to give a particular solution for the potential. 



V{v) = ^FG„(,)e'««r, p(r) = ^/,G„(^)e'«"'-; (A3) y^^^^z) = - P Giz, z')pa,iz') dz' , (A7) 

Gil Gil ^0 Jzo 



thus, the Laplace operator of the Poisson equation 
£oV^V = —p separates into a z and Gm term. 



^^0 



dz^ 



(A4) 



which leads to two equati ons. For Gy =0 (uniform 
field in the vacuum), Eq. (A4) can be solved by inte- 
gration. Using the boundary condition that the poten- 
tial smoothly vanishes at infinity, lim^^oo ^(^) = and 
limz'^oo dzy{z)\z' = the potential in the vacuum is 
given by 



which fulfills the boundary condition that the potential 
vanishes for zq ^ 00 smoothly at infinity. 

If the system has no mirror symmetry, e.g. because the 
field on the top is different from the one at the bottom, 
the charge density can be written a^^^"^ 



'po{z)+p+Y,PGe'^"', \z\<zi 



G^O 



(A8) 



0, 



\z\ > zi 
\z\ < Zi 



P''\^)=\p,{z)+Y.Pae^''-\ \z\>z, (A9) 

Gt^O 
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where —p denotes the average pseudo charge density of 
the interstitial, i.e. 



1 

^^1 J-z^ 



-po,o - ^ po,G^jo{G±zi), (AlO) 



with 



^G„(±N) = -;7rTF^[e^"'^'«G(TNI) 
2eo|(j|||L 

+e-l««l-l^l/?G„(Tkl) 

+e-l««l-l-laG,|(±k|) 



(A14) 



with pG||=o(^) = Z^G^ Po,G^ exjp{iG±z) and Bessel func- 
tion jo{z) = sin z/z. To obtain the potential, one sums 
up the surface charge densities starting from minus infin- 
ity; the electric field is included as surface charge density 
(Jef- For the vacuum region, the Gy = component of 
the potential is given by 



V{z < -zi) 



V{z > zi) 



1 

^0 



a{z') dz' 

) 

[a{z') + (-p)L)vac + ^(-^i)] dz' 



1 

-0{z-zi-^^) 



4f 



(t{z < -zi) = / p{z')dz' 

J — oo 



a{z > zi 



)= /%(^') 

J Z\ 



dz 



(All) 



(A12) 



where V^-^ is the potential at z\ and (/> is a phase due to 
the dipole moment, which is only present if the system 
has neither mirror nor inversion symmetry. The potential 
is given by 



= Vi-z^) + ^ (a(-zl)Z)vac + ^DL(-p)) 



^0 

2i 



y\ ji{G±zi)po^G^zl, 



(A13) 



where ji{z) = (sinz — z cos z)/z^. 

For the Gy 7^ component the solution is given by 
Eq. ( |A7| ); for the FLAPW basis, the integration has to 
be split into two vacua and the film region; for a given 
z > \zi\ the integral in the vacuum can be further split 
into J^^ and J^^"^, which has been done to obtain the 
functions a and P below; the vacuum potential is then 
given by 



n±oo 

aG„(±|^|)=± / PG„(^')e-l°"l-l^'ldz' (A15) 
/3g„ (±1^1) = ± r PG, {zV' dz', (A16) 

J ±Zi 

where pG,, {z) = pa^^iz) ^ S{\z\ - ^^^j^Tcij and (t^^^'^gI 
is the (inhomogeneous) surface charge density of, respec- 
tively, the top and bottom charge sheet. 

2. Dirichlet boundary conditions 

Contrary to the Neumann boundary condition used 
above, which defines the boundary via a surface charge 
density, the Dirichlet boundary condition has a fixed po- 
tential at the boundaries, which matches metallic plates 
held a certain voltage. The Dirichlet boundary condi- 
tions not only imply a different single-particle poten- 
tial but due to the density-dependent image charges 
the effec tive C oulomb interaction between electrons is 
modified."^^'^^" The effect of the latter is not included in 
the described scheme and would require a modified func- 
tional. 

For the constant-potential boundary condition, the 
G|| =0 part, solved by integrating Eq. (A4), is given 

by 



Vo{z) = -- f f po{z")dz"dz' 

^0 J Zq J Zq 

^ — {z - zo)cro + 1^0 
^0 

= -- / a{z')dz' ^-{z-zo)ao^Vo, (A17) 

^0 Jzo ^0 

with a{z) = J^^ p{z') dz' . The lower boundary is set to 
— Zcr, the location of the metallic plate with the associated 
potential Vb = VL^^- For the other plate at +Zcr7 the 
potential Vz^ is given by 



thus. 



V^{z^) 



(JQ, _ I 

60 2zo 



1 

^0 



a(zOdz' + 2z,^ + Fo; (A18) 

^0 



K„ - V- 



- / G{z')dz' 
^0 J-z„ 



(A19) 
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which can be used in a similar way to the Neumann so- 
lution described above. 

To solve the Poisson equation for Gy 7^ 0, we 
start with the solution of the homogeneous problem 
[dl - G||^] Vg|| {z) = 0, which is given by 



which vanishes at the boundary ^z^; for z^r ^ 00 it 



simplifies to Eq. (A6). Using this Green's function, the 
particular solution is given by 



Vg„ (^)=Cie-l^iil^ + C2e+l^iil^ 



(A20) 



The coefficients Ci and C2 have to be chosen such 
that the boundary condition Vg^^ (~^cr) = ^-^J' 
VG|i(^cr) = " cire fulfilled; one obtains 



Ci 



,C2 



; (A21) 



hereby x = exp(|G|| |zcr) and y = exp(— |G|| \za) . On the 
other hand, a Green's function solving Eq. (A4) is given 



G{z,z') 



1 



(A22) 



|G|||sinh(G||2^,) 

J sinh G|| {za- — z') sinh Gy {z + Zcr)-, z < z' 
I sinh Gil (z' + Zcj) sinhG||(zcr — z > z' 



-I' 

£0 J 



V^^''\z) = ^ / G{z,z')po,{z')<iz'. (A23) 



Combining the regular solution of the homogeneous sys- 
tem with the particular solution gives 



Vg,{z) = - r pG,{z')G{z,z')dz' 

-fCie-l^iil^ + Cse+l^iil^ (A24) 

with the Ci and C2 as defined above. 
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